A squeezing property and its applications to 
a description of long time behaviour in 
the 3D viscous primitive equations 

Igor Chueshov* 

Department of Mechanics and Mathematics, 
Kharkov National University, 
Kharkov, 61077, Ukraine 

November 20, 2012 



Abstract 

We consider the 3D viscous primitive equations with periodic bound- 
ary conditions. These equations arise in the study of ocean dynamics 
and generate a dynamical system in a Sobolev H 1 type space. Our 
main result establishes the so-called squeezing property in the La- 
dyzhenskaya form for this system. As a consequence of this property 
we prove (i) the finiteness of the fractal dimension of the correspond- 
ing global attractor, (ii) the existence of finite number of determining 
modes, and (iii) ergodicity of a related random kick model. All these 
results provide a new information concerning long time dynamics of 
oceanic motions. 

Keywords: 3D viscous primitive equations, periodic boundary con- 
ditions, global attractor, finite dimension, random kicks, ergodicity. 
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1 Introduction 



We deal with the 3D viscous primitive equations which arise in geophysi- 
cal fluid dynamics for modeling large scale phenomena in oceanic motions. 
These equations are based on the so-called hydrostatic approximation of 
the 3D Navier-Stokes equations for velocity field u which also contains a 
rotational (Coriolis) force and coupled to thermo- and salinity diffusion- 
transport equations (see, e.g., the survey [21] and the references therein). In 
this model a small variation of the density p of the fluid is taken into account 
via the buoyancy term only and has the form p = po — aT + /3s, where T 
is the temperature and s denotes the salinity. If we take into account one 
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of the factors (T or s) or the diffusivities for T and s are equal, then their 
effect in the dynamics can be represented simply by the density p, or equiv- 
alently the buoyancy b = po — p = aT — /3s (see [23] ) . This is why we prefer 
to deal with the variables u and b. In order to simplify our mathematical 
presentation we impose on u and b the periodic boundary conditions of the 
same type as in [23] (see also [23] and the references therein). Thus we 
consider our problem in the domain 



and denote the spatial variable in O by x = (x,z) = (x\,X2,z) £ O. We 
use the notation V, div and A for the gradient, divergence and Laplace 
operators in two dimensional variable x = (x%, x-z). Below we also denote by 
A XjZ the 3D Laplace operator. The notations S7 x ,z and div XjZ have a similar 
meaning. 

We consider the following equations for the fluid velocity field 

u = (v(x, t); w(x, t)) = (v (x, t);v 2 (x, t); w(x, t)), x = (x,z), 
for the pressure P = P(x, t) and for the buoyancy b = b(x, t): 



v t + (v,V)v + wd z v-vAv-vd zz v + fv ± + VP = G f in Ox(0,+oo), (1) 



b t + (v, V)6 + wd g b - uAb - ud zz b = G b in Ox(0,+oo), (4) 



where v > is the dynamical viscosity, Gf and Gb are volume sources, and 
/ is the Coriolis parameter. We denote v 1 - = (— v 2 ^ 1 ). We supplement 
(fTJ~(j3]) with the periodic boundary conditions imposed on (v;w), P and b. 
We also assume v and P are even with respect to z and w and b are odd 
with respect to z (hence w|^ =0 = and b\ z=0 = 0). These requirements and 
also ([2]), ([3]) lead to the following relations 



w(x,t) = - divv(x, €,t)cl£, P{x,t) = p(x,t) + b(x,£,t)d£ (5) 



for every x = (x, z) € O. It is important (see [3] for the basic discussion) 
that the pressure p is independent of the vertical variable z and depends 
on 2D (horizontal) variable x only. Exploring intensively this observation 
the authors of the paper [3] have implemented a new effective approach 
for proving of the global well-posedness for problems like ([I])-®, see the 
discussion below. 

The system in ([I])-(j3]) was studied by many authors for the different 
types of boundary conditions (see the survey |24j and the references therein). 
The existence of week solutions [22] and (global) well-posedness of strong 
solutions were established (see [3] and also [El [TTl El] ) . The existence of 



O = (0,Lx) x (0,L 2 ) x (-La/2, La/2) C M 3 



div v + d z w = in O x (0, +oo) 
d z P = b in 0x(O,+oo), 



(2) 
(3) 
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a global attractor for the viscous 3D primitive equations was proved in [T5] 
(see also the paper [23] devoted to the periodic case). However, to the best 
of our knowledge, the question concerning dimension of this attractor and 
the related question on determining modes (functionals) are still open. We 
note that the author of [15] claimed the finiteness of the dimension in 2007 
with the reference to a forthcoming paper which is not published yet. 

It should be also noted that stochastic perturbations of ([I])-® were 
studied in [9j [131 E] ( see a l so the references in these publications) and the 
paper [10] deals with a random kick forcing of the primitive equations. 

Our goal is to show that the system in (HJ) — d3J) possesses an additional 
regularity and satisfies a squeezing property in the Ladyzhenskaya form (see 
[20} I21j). This property implies several important facts about dynamics of 
the system and, in particular, the finiteness of the fractal dimension of the 
attractor (see Corollary 13. 6p and the existence of finite number of deter- 
mining modes or functionals (Corollary 13. 7|) . The squeezing property also 
plays an important role in the study of turbulent behaviour in some classes 
of dynamical systems with random kicks arising in the fluid dynamics (see, 
e.g., [IHl IH] ) ■ In this paper we apply the squeezing property established and 
also the theory developed in [19] to present a result on ergodicity of ([I])-© 
forced by random bounded kicks with appropriately chosen frequency (see 
Corollary 

It is commonly recognized (see (3J [HI [III [23J [24] ) that the main difficulty 
arising in the study of the primitive equations in dlj) — (SJ) is related to the 
hydrodynamical part. The calculations which should be made to obtain 
appropriate bounds for a buoyancy type variable are either standard or 
repeat the corresponding argument for the fluid variables in the simplified 
form. This is why below we state our results for the full problem (pQ)-([4]), but 
in the proof we concentrate on the fluid equations only. Moreover, for the 
transparency of our argument we even assume that the buoyancy equation 
is absent. This corresponds to the case of the zero buoyancy (b(x,t) = 0) 
which is possible when Gb = and b(x, 0) = 0. 

The paper is organized as follows. 

Section [2] contains some preliminary material. Here we represent the 
primitive equations in some standard form, introduce Sobolev type spaces 
and state the well-known result on well-posedness of strong solutions. In 
Section [3] we state our main results. The key outcome of Theorem 13.11 is 
the existence of an absorbing ball lying in an H 2 type space. This allows us 
not only guarantee the existence of a global attractor but also provides an 
important step in the proof of the squeezing in Theorem 13.51 In this section 
we also apply our main results to prove the finiteness of the dimension of the 
global attractor, the existence of finite number of determining functionals, 
and ergodicity of a random kicks forced model generated by (P)-©. In 
Section [4] we demonstrate the main steps of the proofs of Theorems 13.11 and 
13.51 in the case when the buoyancy variable b is absent. 
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2 Preliminaries 



In this section we rewrite system (pQH(jl]) in the canonical form, introduce 
appropriate Sobolev type spaces and quote the result on the strong well- 
posedness. 

Using the representations in ([5]) we arrive at the system 



v t + (w, V)w 



div vd£ 



d z v -u[Av + d zz v] + fv~ 



-V 



p(x,t)+ / bd£ 



+ G f in Ox(0,+oo), (6) 



h + (v,V)b- 



div vd£ 



d z b-u[Ab + d zz b]=G b in O x (0, +oo), (7) 



supplied with the conditions: 

"L a /2 



div 
and 



L3/2 



vdz = 0; v is periodic in x and even in z, 



b is periodic in x and odd in z. 
We also need to add intial data for v and b: 

v(0) = wo, 6(0) = 6o- 



vdx = 0; (8) 

(9) 
(10) 



We denote by Hp er ((D) the Sobolev space of order s consisting of periodic 
functions such that f fdx = 0. To describe fluid velocity fields we introduce 
the following spaces with s = 0, 1, 2: 



V s = \v = iy 1 -^ 2 ) 6 [Hp er (0)} 2 : v { is even in z, div J " vdz = o| . 

We also denote H = Vq and equip H with L2 type norm || • || and denote 
by (•, ■) the corresponding inner product. The space V\ and V2 are endowed 
with the norms || • ||vi = ||Va; j2 • || and || • ||y 2 = ||A XjZ • ||. 
To describe the buoyancy we need the spaces 

E s = 1 6 G Hp er (0) : b is odd in z} , s = 0, 1, 2, 

equipped with the standard norms. We also denote W s = V s X E s with the 
standard (Euclidean) product norms. 

As it was already mentioned, starting with [3] the global well-posedness 
of the equations in dU and ([7]) was studied by many authors [El [111 [23l El] • 
In this section we quote the result from |23] (see also [23]) on well-posedness 
of strong solutions under the periodic boundary conditions. 
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Proposition 2.1 f[2~5l [21]) Let G f G V and G b G J5 . T/ien /or every 
= (vo; bo) G Wi problem [6\)- [T0\) has a unique strong solution (u(t); b(t)): 

U(t; Uq) = (u(t); b(t)) G C(R + ; W X ) n L 3 (0, T; W 2 ), VT > 0. 

Moreover, this solution continuously depends on time t and initial data Uq 
and generates a dynamical system (St, W±) with the phase space W\ and the 
evolution operator St defined by the relation StUo = U(t; Uo).This evolution 
operator St possesses the following Lipschitz property 



for every T > and U, U* G B(R) = {U : \\U\\ Wj < R}. 

Remark 2.2 Formally [231 [24] do not contain the proofs of continuity of 
StU with respect to t and the Lipschitz property in W\. However almost the 
same calculations as in |15] alow us to show these properties. We also note 
that in the case considered unique global solvability can be also established 
in smoother classes. For instance, if Gf G V m -\, Gf, G £7 m -l and Uq = 
(vq; bo) G W m for some m > 2 then (see [23l [24]) the solution U(t) belongs 
to the class C(M + ; W m ) n ^2(0, T; W m+ i) for every T > 0. This observation 
is important in our further considerations due to possibility to use smooth 
approximations of solutions in our calculations. 

We conclude this preliminary section with the statement of the following 
well-known uniform Gronwall lemma. 

Lemma 2.3 ([25]) Let g, h, y be nonnegative locally integrable functions 
on [to, +00) such that 



3 Results 

We start with the following assertion which partially follows from Proposi- 
tion 12.11 and based on the calculations given in [23] . 



S t U - S t U*\\ Wl < C t ,r\\U - U*\\ Wl , t G [0,T] 




and 




for any t > to, where a{ > are constant. Then 



y(t + 1) < (03 + a 2 ) e ai for any t > t . 
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Theorem 3.1 (Smooth Absorbing Ball) Let Gf G V\, G b G E\ and 

also (d z Gf,d z Gb) € [Lg(O)] 3 . Then the dynamical system (S t ,W\) gener- 
ated by problem Kb])~ /I10\) is dissipative and compact Moreover, there exists 
> such that the ball 

m = {U eW 2 : \\U\\ W2 < R*} 

is absorbing, i.e., for any bounded set B in W\ there is ts such that 

S t B C SB for allt>t B - 

In this ball SB there is a forward invariant absorbing set Qi . 

We note that the dissipativity and compactness of (St, W\) stated in Theo- 
rem [37T] is known from |23j under slightly weaker conditions on Gf and Gb- 
However the existence of an absorbing compact set which is bounded in W2 
requires additional hypotheses and additional calculations. 

Remark 3.2 In the case when Gf = and Gb = one can prove (see some 
argument below in Section [4. ip that 

Ve>0, 3T(R,e): \\S t U\\ Wl < e Vt > T(R,e), \\U\\ Wl < R. 

This means that in the case of vanishing forces the zero equilibrium state is 
asymptotically stable. 

Using the standard results (see, e.g., one of the monographs [TJ G2 [25] ) 
on the existence of global attractors we can derive from Theorem 13.11 the 
following assertion. 

Corollary 3.3 (Global Attractor) Let the hypotheses of Theorem \3.1\ be 
in force. Then the dynamical system (St,W\) generated by problem [h^)~ (10\) 
possesses a compact global attractor 21 which is a bounded set in W^- 

Remark 3.4 The result on the existence of a global attractor for 3D vis- 
cous primitive equations was known before, see [15], and also Remark 2.1 in 
|23j for the periodic case. Our improvement is that we state the bounded- 
ness of 21 in H 2 type Sobolev space. One can also show that in the case of 
sufficiently smooth sources Gf and Gb the attractor 21 possesses additional 
spatial smoothness. Using the method presented in |23] we can even prove 
Gevrey regularity of the elements from the attractor provided the sources 
possess appropriate smoothness properties. However we do not pursue these 
improvements because our main goal is a squeezing property with applica- 
tion to dimension and ergodicity. 

1 For the definitions see [I], for instance. 
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In the space W\ we consider the bilinear form 

a(U, £/*) = / [V X)Z v ■ V x , z v* + V x>z b ■ V^A] dxdz, 
Jo 

where U = (v;b) and C7* = (u*;6*) are elements from W\. This form gen- 
erates a positive self-adjoint operator A with a discrete spectrum, Due to 
periodicity of the problem the corresponding eigenfunctions and eigenvalues 
can be easily described. We assume that they are reorganized in such way 
that 

Aek = Afcefc, < Ai < A 2 < . . . , lim A& = +00. 

Moreover, {e^} is the orthonormal basis in Wq. We denote by P/v the 
orthoprojecton on the Spanjei, . . . , e./v} and Qn = I — Pn- 
Now we can state our main result. 

Theorem 3.5 (Squeezing) Let the hypotheses of Theorem \3.1\ be in force. 
Then the dynamical system (St,W\) generated by problem $b))-[T(J\) on the 
forward invariant absorbing set Qi possesses the properties: 

• Lipschitz property: For any U,U* E we have the relation 

\\S t U - S t U*\\ Wl <C & e as,t \\U -U4 Wl , Vi>0, (11) 
where C@ and a@ are positive constants. 

• Squeezing property: For every T > and < q < 1 there exists 
N = N(T, q) such that 

\\Qn[S t U - S T U*]\\ Wl <q\\U-U4 Wl (12) 

for any U and U* from the absorbing forward invariant set Ql . 

The squeezing property stated in Theorem 13.51 allow us to apply Ladyzhen- 
skaya Theorem on the dimension of invariant sets (see [20} [21] and also 
Theorem 8.1 in [5l Chap.l]). This yields the following assertion. 

Corollary 3.6 (Finite Dimension) Under the hypotheses of Theorem \3.1\ 

the global attractor 21 of the dynamical system (St, W\) generated by problem 
H)\)~ il0\) has a finite fractal dimension. 

This means that the asymptotic dynamics in this model is finite-dimensional 
and topologically can be represented by a compact set in M. d with appropriate 
d < 00. 

The next outcome of the squeezing property is the existence of a finite 
number of (asymptotically) determining modes. We note that this notion 
was introduced in [12] for 2D Navier-Stokes system and was studied by 
many authors (see the discussion in the recent paper |llj). For a general 
theory of the determining functionals we refer to [1] , see also [5] and [7] for a 
development of this theory based on the notion of the completeness defect. 



7 



Corollary 3.7 (Determining Modes) Let the hypotheses of Theorem \3.1\ 
be in force and N be such that f)12|) holds for some T > and < q < 1. 
Then the dynamical system (St,Wi) possesses N determining modes. This 
means that the relation 

lim \\P N [S t U - S t UM Wl = for some U, U* G Wi (13) 

t— s-+oo 

implies that \\StU — StU*\\w! — >■ as i — >• +oo. 

Proof. The relation in (j!2j) implies that 

||S , T Z7- < ST^||m < q\\U-U*\\ Wl + \\P N [S T U-STU*}\\ Wl , 0<q<l, (14) 

for any U and C7* from the absorbing forward invariant set Iterating the 
relation in (j!4j) we obtain that 

n-1 

ll^-S^ET.Hwi <<7 n ||C/-C/*||^+^g n ^ fc ||Piv[^-'5r^]llw 1 (15) 

fc=0 

for every n = 1,2, . . ., where {/ and C/* are from '3. This inequality and 
also the Lipschitz property in (fTT|) makes it possible to prove that relation 
(fT3]l implies that \\StU — StU*\\wi — > as t — > +oo. For some details in the 
abstract situation we refer to the proof of Theorem 1.3 in Chap. 5]. □ 

Remark 3.8 Using the same idea as in [7] it is also possible to derive from 
(|15p the existence of other finite families of determining functionals. 

Now we consider an example of an application of the results above to 
ergodicity of a random kick forced model generated by problem ()10|) . 
We deal with the simplest situation and suppose that Gj = and Gj, = 0. 
Moreover, we consider a model in some fixed ball from W\ and assume that 
the frequency of kicks is smaller than a certain critical value (depending on 
the radius of the ball and the amplitude of the kicks). We plan to remove 
all these restrictions and discuss further developments in our joint studies 
with S.Kuksin and A.Shirikyan. 

Let {ijk} be i.i.d. random variables in W\. Assume that the support of 
the distribution V(r]i) in W\ contains the origin and is bounded, i.e., there 
exists Rkick > such that 

suppX>(r/i) C B{R kick ) = {U G W\ : \\U\\ Wl < R kick } ■ 

Here and below we keep the notation B(R) for the ball with the center at 
of the radius R in the space W\ . 

Let us fix R > R k i ck and choose T c = T(R, R k i c k> ^) > 1 such that 

S t B(R) C B(R - R kick ) and S t -xB{R) C 9 for all t > T c , 
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where Q is the forward invariant absorbing set given by Theorem 13. 1[ This 
choice of T c is possible due to Remark 13.21 and Theorem 13.11 

Now we take T >T C and define a Markov chain in W\ by the formula 

U k = F k {U k ^ 1 ) = S T U k ^ 1 + r lk , k = 1,2,... (16) 

This chain is called a random kicks model generated by ([6])- (|10p (with Gf = 
and Gb = 0). The parameter T _1 has meaning of the frequency of the 
kicks. We refer to [19] for motivation an physical importance of different 
types of kick models in the turbulence theory. 

It is obvious that if Uq 6 B(R), then we also have that U k £ B(R) for 
all k. Thus we also have a chain in B(R). 

It follows from Proposition 12. II that 

\\St-iU - St-iU*\\wi < C t ,rW U ~ U *\\W!, for all U, U* € B(R). 

Let us fix < q < 1 such that r\ = qC T ^ < 1. Since St-\U, St-iU* S 2, 
we can apply Theorem 13.51 and choose N such that 

\\Qn[StU - S T U*\ \\ Wl < q\\ST-iU - S T -iU*\\ Wl < v\\ u - U*\\ Wl 

for all U, f7* € B{R). Thus we have the squeezing property on the ball B{R). 

Applying now Theorem 3.2.5[19j and the properties of the mapping St 
established above we arrive to the following assertion. 

Corollary 3.9 (Ergodicity) Assume that Pntji and Qn^i o.re indepen- 
dent and the distribution of P^Tji in the (finite- dimensional) space PnW\ 
has a density p(U) with respect to the Lebesgue measure dU on PnW\ such 
that 

[ \ P (U + V)-p(U)\dU <C\\V\\ Wl , W€P N Wi, 

JP N W! 

where C is a constant. Then for the chain H6\) there is a unique invariant 
probability Borel measure fi on B(R). Moreover, there is < 7 < 1 such 
that 

||F fc *A-/42<C0R,A)- 7 fc , k = l,2,... 

for any probability Borel measure A on B(R), where \\ ■ \\* L denotes the cor- 
responding dual Lipschitz norm (see, e.g., 

To obtain this result we apply |19|. Theorem 3.2.5]. Conditions 3.2.1- 
3.2.3 in [19] follow from the properties of St stated above. We also note that 
[19^ Theorem 3.2.5] formally deals with a Markov chain in a whole space. 
However one can see that the arguments from |19| remain true if we restrict 
the iteration procedure in (|16p on an invariant subset (this observation was 
already used in [6] to prove ergodicity for random kick-forced 3D Navier- 
Stokes equations in a thin domain). Another possibility to fix this issue is 
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to consider the extension St of St from B(R) on the whole space W\ by 
the formula 

' S T U, \\U\\ Wl < R; 

S T U = < , ~ N 

I \\u\\ Wi >r. 

One can see that St maps W\ into B{R — Rkick)- Thus after the first step 
in ()16p we are in B(R) and therefore we have St = St for k > 2 in (|16p . 
We also note that St is globally Lipschitz (the proof of the latter property 
can be found in [8j p. 64], for instance) and the Conditions 3.2.1-3.2.3 [19] 
follow directly. 



4 Proofs 

Our arguments are more or less standard and use the methods developed in 
[31 [151 123] (see also the survey [21] ) . In fact they are some refinement of the 
calculations known now from [3] , see also |15[ \23\ [24] ) . This why we consider 
only key steps in the argument. Moreover, we deal with the reduced system 
which appears in the following way. 

If we assume Gb = and bo = 0, we can take b = as a solution to ( [ 7 ] ). 
Thus we arrive to the problem: to find a (horizontal) fluid velocity field 

v(x,t) = (v (x, t); v 2 (x, t)), x = (x,z), 

and the pressure p(x, t) satisfying the equation 



v t + (v, V)u 



div vd£ 



d z v - v[Av + d zz v] + fv^ = -Vp + G f (17) 

in O x (0, +oo) supplied with conditions ([8]) and with initial data 

v(x, z,0) = vo(x, z). (18) 

It follows from Proposition 12.11 that for every Gf G V\ and v € V\ problem 
(fT7|) , (fT8j) , (j8]) has a unique strong solution 

v(t) eC{R + ;V 1 )nL 2 (0,T;V 2 ), VT > 0. 

These solutions generate a dynamical system (St, V\) in Vi. Below we prove 
our main result for this system. Thus we deal with calculations involving 
the velocity field v, but not the buoyancy. Moreover using approximation 
procedure and the results from |23[ [24"] we can even assume that all our 
calculations deal with smooth solutions. 



4.1 Preliminary a priori estimates 

In order to prove Theorem 13.11 for the system (St,Vi) generated by (fT7|) we 
need only to check the existence of absorbing ball in the space V%. For this 
we need to repeat and partially refine some calculations in |15[ [23] . 
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4.1.1 A priori estimate in H 

If we multiply (fT7|) by v in H, then we obtain 

d|M 12 



dt 

which implies that 

rt+l 

\v> + "' 2 



+ v[\\Vv\\ 2 + \\d z v\\ 2 ]<c\\G f \\< 



+ u J [||V«|| 2 + ||<9^|| 2 ]dr < 2e- cot ||i»(0)|| 2 + c 1 ||G / || 2 

for all t > 0. Thus there exists Rq > such that for any i? > there is 
tR> such that 

IK*)H 2 + ^ [l|V«|| 2 + ||5^|| 2 ]dr < Rl for all i > i fl (19) 

with ||f (0)|| < i?. Moreover, we can choose i? 2 , = e + ci||G/|| 2 with arbitrary 
e > 0, in this case i# also depends on e. 



4.1.2 Splitting 

Let 



v = — vdz = (v) z and v = v — v. 

^3 J-L a /2 

As in [3] one can see that the fields v and v satisfy the equations 

v t + (v, V)v + (v, V)v + (V, v)v - uAv + /u" 1 = -Vp + ( 20 ) 
with div u = in T 2 = (0, Li) x (0, L 2 ), and 



z 

div vdz 



o 



d z v 

+ (£ , V)u + (v, V)v - [(v, V)v + (V, v)v] 

-v[Av + d zz v] + fv x = G f (21) 

inO = T 2 x (-L 3 /2,L 3 /2). 
4.1.3 .ff^-estimates 

Now we multiply (|2ip by |v| 4 i) and integrate over 0. In the same way as in 
[31 [15] we obtain 



d "~" R + v J [\Vv\ 2 + \d z v\ 2 ] \v\ 4 dx 

< C [\\v\\ 2 \\Vv\\ 2 + \\Vvf] \\v\\l 6 + C.WGfUMW,- 
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This implies that 



d 



J t Ml 6 < Co [1 + IMI 2 ] IIIIVHI 2 • NIL + ClHG/lklHk 
< Co [l + (l + lbll 2 )lll|VH| 2 ]-||^lli 6 +C71IIG/III,. 



Since ||£||l 6 < c[||V?;|j 2 + ||5 z f;|| 2 ], we can apply Lemma 12.31 and relation 
(|19p to obtain that there exists R\ > Rq > such that for any R > there 
is t* R > tR such that 

\\v(t)\\ L6 < R x for all t>t* R with \\v(0)\\ < R. (22) 

Moreover using (HI]) with Rl = E + c 1 \\G f \\ 2 we can conclude that 

Ve€(0,1]: \\v(t)\\u <C[e+\\G f \\ 2 ] for all t > t* R>£ 

with ||v(0)|| < R. Using (|22p we can also assume that 
ft+i f 

/ / [\Vv\ 2 + \d z v\ 2 ] \v\ A dxdr < C(Rt) for all t > t* R (23) 
Jt Jo 

with ||f(0)|| < R. In the case Gf = we can change C(Ri) into Ce and t R 

Next we multiply (j20|) by — A-u in L2(T 2 ). As in [3j using the relations 

/ (v, V)vAvdx = 0, / VpAvdx = 0, / v ± Avdx = 0, 
Jt 2 Jt 2 Jt 2 



we have that 



(v, v)*j + (y,v)v 



Avdx 



J2 



GfAvdx. 



This implies (see [3] for some details) that 



^\\Vv\\ 2 + v\\Av\\ 2 <C\\\Vv\\ 2 + I \v\ i \Vv\ 2 dx + \\Gf\\' 



o 



Therefore ([22]) and (|23j) give us that there exists R2 > such that for any 
-R > there is t* R * > t* R such that 

rt+l 

\\Vv(t)\\ 2 + u J \\AvfdT <R 2 2 for all t > t* R * with ||u(0)||<JJ. (24) 

As above in the case Gf = we can take R\ of the order e (with t R 
depending on e) . 

The next step is the estimate for v z . We first note that u = v z solves 
the problem 



ut + (v,V)u 



div vd£ 



d z u - u[Au + d zz u] 
= - fu L - (u,V)v + (V,v)u + d z G f (25) 
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in the class of periodic (odd in z) functions. Since u u = and 



(v,V)u 



div vd£ 



d z u > udx = 0, 



using the multiplier u after integration by parts we obtain that 
ld\\u\\ 2 



2 dt 



+ ^[||Vn|| 2 + ||<9 2 u|| 2 ] < ||5,G/|||H| + / M|u||V 

Jo 



u\dx, 



which implies (for some details we refer to [3]) that 
d\\uf 



dt 



+ ^[||Vu|| 2 + ||cU|| 2 ] < C [IMUJMI 2 + \\d z G f \ 



<C[ll^lli 6 + l|V^|| 4 J \\u\\ z + C\\8 z G f \ 



(3) 

Therefore there exists R3 > such that for any R > there is > t R 
such that 

ft+l 

+ v J [\\Vv z f + ||^|| 2 ] dr < Rl for all t > vjjj,' 

with ||f (0)|| < R. Moreover, R 2 ~ e in the case Gf = 0. 

Next we estimate ||Vt> ||. For this we multiply (fT7]) by — Av. As in [3] we 
obtain 

d\\Vv\\ 2 



dt 



+ v[\\Av\\ 2 + \\V8 Z 



<Cni^l|l 6 + llV7;|| 2 ||^|| 2 l ||V V || 2 + C||G f 112 



which implies again that there exists R4 > such that for any R > there 
is tft > tfi* such that 



\Vv 



for all t > t$ with ||v(0)|| < R. Moreover one can see from the analysis 
above that in the case Gf = for any < e < 1 there exists i_R j£ > such 
that 



+ v 



t+i 



\Av\\ 2 + ||V^|| 2 + H^ll 2 ] dr < R\ (26) 



\Vv 



rt+l 

+u J [\\Av\\ 2 + \\Vv z \\ 2 + \\v zz \\ 2 ] dr <e for all t > t R>E (27) 



provided ||f(0)|| < R. 

The relation in (j26j) yields that the system (St, V\) is dissipative in V\. In 
the case Gf = by (127p an absorbing ball can be chosen with an arbitrary 
small radius which implies that the zero equilibrium is asymptotically stable 
in this case (as it is claimed in Remark 13. 21 when the buoyancy is presented). 
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4.2 Completion of the proof of Theorem 13.11 

To conclude the proof we need to show uniform smoothing of trajectories. 
We note that the result for individual trajectories is known from [23] (see 
also [24] )• However in these references the size of the corresponding ball is 
not controlled and thus we need some refinement of that results. This is 
why below we mainly follow the line of argument in [23] (see also [24J ) . 
We multiply (HZD by (-A x ^ z ) 2 v to get 

- — \\A x , z v\\ 2 + u\\[-A x ^ 2 v\\ 2 = - y i( v > v ) v + w(v)d z v] {-A x , z ) 2 vdx 

+ (G f ,(-A XtZ ) 2 v), (28) 
where w(v) = w(x,t) is defined by ([5]). We also use the facts that 

/ v ± (-A XyZ ) 2 vdx = I [A XjZ v] ± A x , z vdx = 
Jo Jo 

and 

/ Vp(-A x , z ) 2 vdx = 0. 
Jo 

The general structure of the nonlinear terms in (|28p is given by the formula 

[(v, V)u + w(v)d z v] (-A XjZ ) 2 vdx 

A x , z [(v, V)v + w(v)d z v] A x>z vdx 

■■ I [(v, V)A x , z v + w{v)d z A x>z v] A x>z vdx + E(u) = E(«). 
Jo 



N(v) 



Here 



dx, 



(29) 



S(w) - V / D 2 v i D 1 v j D 2 v j + D k w(v)D 2 - k d z v j A x . 

where k,i, j = 1,2, D k is a differential operator of order k, and the sign ~ 
means that we deal with linear combination of the terms in the right hand 
side. 

Using the embeddings H x (0) C L 6 (0) and H l / 2 {0) C L 3 (C7) (see [26]) 
and also interpolation one can see that 



Si = 



[ D 2 v t D 1 v j A x , z v j dx 
Jo 



< C||D 2 ^IH| j D 1 ^||l 6 || j D 2 ^ 



< C||^|[l, 2 |l J D 2 ^|| 1 /2 | | jD2 ^||V^ < e ||[_ Aa; , ] 3 / 2 T;|| 2 + 
To estimate the second term in 



1 10/3 



we start with the case k = 2: 



D 2 w(v)d z v j A x ^dx < C\\D 2 w(v)\\\\d z vi\\ L jD 2 v j \\ L3 

< C \\ v \\h3 \\dzV\\L 6 \\v\\H2 

<e\\[-A x , z f/ 2 v\\ 2 + C £ \\d z v\\ijv\\ 2 H2 . 
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To estimate the term with k = 1 in (|'29p we use the following lemma (see 
Proposition 2.2 in [2]). 

Lemma 4.1 For smooth functions (f> and tp and vector field v = (v l ]v 2 ) we 
have 




divv(x, £,)dl; 



ifxfidx 



<c\\v\\fAvCAnfM ll2 M 



This lemma yields the following estimate 



f D l w{v)D l d z v^ x ^dx < C||^||^ 2 |kll^ 2 ||^^ / 2 2 ||^'||J/ X 2 ||L> 2 ? ;| 
Jo ... 



<C\\v\\ H 3\\vf H2 <e\\[-K x>z f 2 v\\ 2 + C £ \\v\\m- 
Thus we arrive at the following relation 



d 



-IIA^H 2 + v\\[-A x , z ] 3 / 2 v\\ 2 <d [1 + ||A^|| 2 + II^IIIJ ||A^| 



+ C2\\G f \\ 



(30) 



Our next step is an estimate for ||<9 2 t> ||l 6 . As above let u = d z v. We multiply 
([25]) by \u\ 4 u. Using the relations u^u = and 



o 



{v,V)u 



div vd^ 



d z u > \u\ A udx = 0, 



we obtain that 



1 d 1 1 ,i6 
6dt U 111/6 



o 



[\Vu\ 2 + \d z u\ 2 ]\u\ 4 dx 



\u\ u). 

Jo 



In the same way as in [231 [24] after integration by parts we have that 

/ {-(u, V)v + (V,v)u} \u\ 4 udx < C\\v\\ Lao / \Vu\\u\ 5 dx 
Jo Jo 



< C\\v\ 



Jo 



u| 2 |u| 4 <ix 



1/2 



U/ll 3 
\ u Wl 6 



<e [ \Vu\ 2 \u\ 4 dx + C £ \\v\\ 2 L J\u\\ 6 Lti . 
Jo 



This implies that 



4lHli 6 +^ [ [\Vu\ 2 + \d z u\ 2 ]\u\ A dx 

at Jo 

<c[||«||| r ,|Hli 8 + ||a*G / || Lfl Hi 8 ]. 
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Thus we have that 



d 

~dl [ 



u 



<d [l + \\A Xr 



u 



li 6 + C 2 \\d z Gf\\ L6 . 



Using Lemma 12.31 and the estimate in (|26p we obtain that 

\\d z v(t)\\ Le < R 5 for all t>t$ > t$ with ||u(0)|| < R. 

Now returning to (|30p and applying Lemma [2.3l allow us obtain the following 
assertion. 

Proposition 4.2 Let Gf G V\ and d z Gf € Lq(0). Then there exists R* > 
such that the ball SB = {v £ V 2 : \\v\\v 2 5; R*} is absorbing for the 
dynamical system (Sf,Vx) generated by problem Moreover, there is 

a forward invariant absorbing set & in this ball. 

Thus the proof of Theorem 13. II (in the case bo = and G^ = 0) is complete. 
4.3 Squeezing estimate 

Now we prove Theorem 13.51 (in the case bo = and Gb = 0). We start with 
the following assertion. 

Proposition 4.3 For any v\,V2 from the absorbing forward invariant set 
& we have the relation 

\\Sm- S t v 2 |k <C 9 e aat \\v 1 -v 2 \\vi, Vi>0, 

where C$> and are positive constants. Moreover, for v*(t) = §tv* and 
v**(t) = Sty** we also have that 



K*)llvi+/ \H T )\\v2 dT ^ Coje^Wv, -v*4 2 Vl , Vi>0, 
Jo 

where u(t) = v*(t) — v**(t). 

Proof. We use the same type of argument as in [15]. For this we write the 
equation 



u t - v[Au + d zz u] + (u, V)t>* + V)u 



div ud£ 



d z v* 



div v**d£ 



d z u + fu 1 - = -V(p* - p**) (31) 



for the difference u(t) = — i>**(i) = Stu* — Sty**. Then we multiply 
([3Tj) by (— A X)Z )u and apply the same argument as in [15]. At the final stage 
we use the fact that 



II U *WI|V2 + Ik**(*)l|v2 < Ccj, t > 0, 

for initial data v* and v** from 



□ 
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Let {Afc} be the eigenvalues of the corresponding Stokes type operator 
A which arises in problem (|17p and {e^} be the corresponding eigenvectors. 
We can assume that 



< Ai < A 2 < . . . 



lim Afc = oo. 

k— >oo 



We denote by P/v the orthoprojector in H on Spanjei, 
/ — Ppf. We multiply (|3T|) by QnAu and obtain that 



, e^}. Let Q 



N 



l£\\A^Q N u\\ 2 +v\\AQ N uf 



-([(«, V)w* + (v**,V)u],QnAu) 



+ 



+ 



o 



o 



div ucf£ d z v*QNAudx — f(u ,Q^Au) 
divv**d!; dzuQ^Audx. (32) 



Now we need the following lemma which follows from the Holder inequality 
and from the standard Sobolev embedding (see, e.g., |26j): 

H S (0)CL P (0) when s - 3/2 = -3/p, p>2, OcR 3 . 



Lemma 4.4 Let Oct 3 and u, v be sufficiently smooth functions. Then 
for any < s, r < 3/2. such that r + s = 3/2. 



Using this lemma to estimate the right hand side in (|32p on the forward 
invariant set & we can find that 



1 d_ 

2dt 



^^QNuf+vWAQNuf < C Rt \\u\\ H3/2 \\AQ N u\ 



where i?* is the same as in the statement of Proposition 14.21 This yields 

j t \\A l l 2 Q N uf+v\\AQ N u\\ 2 < C R ^ e \\u\\ 2 Vl +s\\u\\ 2 V2 
for any e > 0. This implies that 

WA^Qnu^W 2 < e- vXN+lt \\A l / 2 Q N u{>d)\\ 2 

+ C Rtt£ [ e-" Ajv + 1 (*- T ) dr + e / e -^w+i(*-r) \\ u \\ 2 V2 dT. 

J J 



Using Proposition 14.31 we obtain that 



{A^Q 



N u 



< 



-i/Ajv+i* _|_ 



Cr, 



OR t + ^AjV+l 



+ eC R , e 



,OlR,t 
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Now for any fixed T > and < q < 1 we can choose e and N such that 

\\A l l % Q N u{T)\\ 2 <q 2 \\A l / 2 u{Q)f. 

Thus we arrive to the squeezing property: For every T > and < q < 1 
there is N such that 

\\Qn[Stvi - S T V2]\\vi < q\\vi - v 2 \\vi 
for any vi,V2 from the absorbing forward invariant set . 
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